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AN INVESTIGATION OF SCAVENGING OF 
RADIOACTIVITY FROM NUCLEAR DEBRIS CLOUDS: 
RESEARCH IN PROGRESS 


Abstract 


The publication Potential Exposures 
from Low-Yield Free Air Bursts, 
addressed the problem of rapid deposition 
of radioactivity from nucleai- debris 
clouds at various distances from ground 
zero, 
The results showed that large 
values of infinite exposure and large 
exposure rates could occur at the surface 
under the following special circumstances: 
the radioactivity is distributed in the 
lower troposphere, it is scavenged at an 
early time by a rain-bearing system and 
rapidly deposited over the ground, and it 
is confined to a relatively small area. 
More recent studies have been and are 
being made at LLL to make it possible 
to assess 1) the reasonability of the r e ­ 
sults of this earlier publication by con­ 
sidering in greater detail the physical 
processes involved and 2) the extent to 
which the results may be changed by man- 
made modifications. This paper is a prog­ 
r e s s report on the research projects cur­ 
rently underway at LLL which are directed 
toward quantifiable assessments for client 
needs. Specifically, these projects are: 


1) A study of the dynamics of cumulus 
cloud interaction with a nuclear 
debris cioud to determine possible 
enhancement or reduction in peak 
surface concentrations. 
2) A study of the microphysical inter­ 
actions of debris particles with 
cloud droplets and raindrops to 
determine rainout and washout 
eff ic iencies. 
3) A study to determine the feasibility 
of rendering the debris particles 
non-wettable, thereby effectively 
making the rainout coefficient 
negligibly small. 


4) A study of the potential effects of 
dry deposition due to the detonation 
of larger numbers of low yield 
devices within a confined area. 
5) The development of a three- 
dimensional atmospheric diffusive 
particle-in-cell code capable of cal­ 
culating the time-dependent distri­ 
bution of air pollutants under con­ 
ditions of turbulent diffusion and 
wind shear. 


Introduction 


A major conclusion from the report 
Potential Exposure from Low-Yield Free 
Air Burst" was that the hazard of close- 
in exposure from these bursts may be 
appreciable when scavenging by precipita­ 


tion occurs. For example, rainout of the 
radioactivity within the debris cloud of a 
1-kt nuclear detonation at 100 km down­ 
stream of ground zero may result in 
external gamma dose rates as high as 
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100 R/hr and infinite whole body doses of 
1000 R or higher. A study by Conner and 
2 
Phillips further substantiates the hazard 
inherent in these low-yield bursts. 
From the interest generated by these 
papers, it has become apparent that a 
greater effort is required to obtain more 
detailed and more accurate assessments 
of the hazards involved. With this need 
in mind, the Atmospheric Sciences Group 
at Lawrence Livermore Laboratory some 
months ago initiated work on a suite of 
research projects to improve understand­ 
ing and assessment capabilities pertain­ 
ing to the scavenging of debris from 
nuclear clouds. It is the intent of this 
paper to present an outline of eaeh of 
these projects, their current status and 
some preliminary results. 


The specific research problems to be 
discussed are as follows: 
(1) A model of the interaction between 
a CUIT'VIUB cloud and a nuclear debris 
cloud to determine possible enhancement 
or reduction in surface concentrations 
due to relative positions and motions. 
Among the assumptions used in the early 
calculations was that the vertical integral 
through cloud-center was deposited on the 
ground and integrated exposures and 
exposure rat ,s were calculated based 
upon this deposition. It is possible, how­ 
ever, that debris will be entrained into 
a developing cumulus system and that the 
resultant vertical integral of the debris 
concentration in the cumulus cloud may 
be enhanced by this redistribution of 
radioactivity. Conversely, the resultant 
vertical integral may be decreased de­ 
pending upon the relative vertical positions 
of the scavenging system and the debris 
cloud. 
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(2) A study of the microphysical inter­ 
actions of debris particles with cloud 
droplets and raindrops to determine 
deposition rates. It is known that the 
rate at which the debris is removed from 
a cloud depends upon the rain rate, the 
spectrum of particle and raindrop sizes, 
the relative fall velocities of the raindrops 
and the particles, the efficiencies of 
capture for a given particle size, and the 
chemical properties of the particles. 
Frcm these parameters washout and 
rainout efficiencies can be calculated 
and some measure of the effective­ 
ness of the rainout process can be 
determined. 
(3) A study to determine the feasibility 
of rendering the debris particles non- 
wettable. Particles resulting from free 
air bursts appear to be insoluble and 
wettable. This being the case, particles 
with radii as small as 0.1 jum can serve 
as condensation nuclei. According to 
3 
arguments by Fletcher and experimental 
results by Twomey, particles whose 
surfaces depart from being completely 
wettable cannot serve as condensation 
nuclei under typical cloud super saturations. 
Since nucleation scavenging appears to be 
the dominant rainout mechanism, a slight 
"eduction in t 
: e particle wettability could 
greatly reduce the potential hazard from 
rainout. 
(4) A study of the potential effects cl 
dry deposition (deposition caused by the 
dispersion of debris particles to the sur­ 
face) due to the detonation of a large 
number of low yield devices within a con­ 
fined area. 
(5) The development of a three- 
dimensional Atmospheric Diffusive 
Particle-in-Cell Code capable of 


calculating the time-depending distri­ 
bution of air pollutants under conditions 
of turbulent diffusion and wild shear. 
A summary of the pertinent findings of 


the report Potential Exposures from Low- 
Yield Free Air Bursts is included for 
the benefit of readers not familiar with 
this previous report. 


stabilization altitude. 
5. The debris cloud undergoes no 
shear. 
6. The debris cloud is dispersed 
radially and vertically by diffusive 
eddies in the atmosphere. 
7. At some point downstream, the 
debris in the debris cloud is sub­ 
jected to scavenging by a rain 
bearing system, thereby causing 
all or a portion of the radioactivity 
to be deposited upon the ground. 
This deposition process is assumed 
to be instantaneous. 
8. Once the radioactivity is deposited 
on the ground, it remains in place. 
Assuming the above sequence of events 
occurs, it is possible to calculate a rate of 
exposure and an infinite exposure to a per­ 
son standing at the depositionpoint once a 
means of determining the areal concen­ 
tration of the radioactivity is specified. 
For this study, the areal concentra­ 
tion is defined as the radioactivity re­ 
moved from a column of unit cross 
sectional area through the center of the 
debris cloud (also often referred to here 
as the vertical integral ). Thus as the 


Depending upon use, the term vertical 
integral may refer to the radioactivity 
deposited, the dose rate due to this radio­ 
activity, or the infinite dose due to this 
radioactivity, and, unless otherwise 
specified, it refers to cloud center. 


Summary of the Pert 
Potential Exposures from 


The risk due to fallout from a free air 
burst has generally been considered 
negligible since few particles large 
onough to settle to the ground are created. 
Early time rainout or washout of radio­ 
active debris from air bursts has been 
considered a potential hazard^ but since 
choice of meteorological conditions has 
always been available during testing, the 
probability of dpmage from this phenomenon 
has been consicered small or non-existent. 
If, however, it becomes necessary to 
detonate many low yield devices within a 
small area, at a time when choice of 
meteorology is not an option, then hazards 
due to rainout or washout become much 
more probable. 


Consider the following assumptions 
relating to the detonation of an all fission 
nuclear device: 
1. It is detonated at a low altitude but 
high enough that there is no inter­ 
action of the fireball with the sur­ 
face (free air burst). 
2. The resulting debris cloud rises to 
its stabilization altitude and is 
approximated in form by a right 
circular cylinder. 
3. The distribution of radioactivity in 
the debris cloud at stabilization is 
assumed Gaussian in both the radial 
and vertical directions. 
4. The debris cloud moves downstream 
at the speed of the wind at the 


Summary of the Pertinent Findings from 
Potential Exposures from Low-Yield Free Air Bursts 


ft. 


debris cloud moves downstream, its 
vertical integral can be calculated as a 
function of time from the diffusion param­ 
eters, and the pertinent exposures and 
rates can then be calculated. The vehicle 
for calculating the vertical integral is the 
Large Cloud Diffusion Code 2BPUFF. 
6 


Table 1 represents the input param­ 
eters to 2BPUFF for vertical integral 
calculations of debris clouds from one, 
ten, and one-hundred kt detonations. The 
cloud dimensions, height, or r i s e , 
7 and 
meteorology 
are based upon mean con­ 
ditions. Two diffusion conditions, 
atmospheric dissipation (e> and vertical 
diffusivity (K ), are assumed in order to 
obtain a reasonable range of values of 


the vertical integral as a function of 
6 9 10 
time. ' ' 
Figure 1 shows the vertical 


integrals {converted to time of arrival 
dose rates) in R/hr for the three yields, 
and Fig. 2 shows the infinite whole body 
exposures. 
The time-dependence has 
been converted to dependence on distance 
from ground zero in both /igures. 
It is 
seen that extremely large exposures and 
rates can occur at distances even as far 
as 1000 to 2000 km for the 100 kt case. 
However, another important factor is that 
the debris clouds for the large-yield 
bursts will rise to greater heights and 
somr or all of the debris will no longer 
be available for scavenging because no 
rain bearing systems will reach high 
enough. Thus, based upon an estimate 
of types of scavenging systems, it has 
been assumed in this report thai 100% of 
the 1-kt debris, the lower 10% of the 


Table 1. Input cloud parameters for exposure calculations. 


Center 
height 
(m) 


Thick­ 
ness 
3 
(m) 
Radius 
3 
to) 


« 
b 
Gross 
fission 
(pCi) 


131, 
(pCi) 


Wind. 
speed 
(km/hr) 


Atmospheric 
dissipation, 
e (ergs/g-sec) 
Vert icaldiifusivity 
(cm 
2/uec) 
yield 
<kt) 


Center 
height 
(m) 


Thick­ 
ness 
3 
(m) 
Radius 
3 
to) 


« 
b 
Gross 
fission 
(pCi) 


131, 
(pCi) 


Wind. 
speed 
(km/hr) 
Slow 
Fast 
diffusion 
diffusion 
Slow 
diffusion 
Fast 
diffusion 


1 
10 
100 


2,840 
7,000 
11,700 


1760 
3060 
5340 


920 
2400 
6000 


4.44 X ID 
2 0 


4.44 X 10 
2 1 


4.44 X 10 
2 2 


1.4X 10 
1 7 


1.4 X 10 
1 8 


1.4 X 1 0 
1 9 


39.6 
70.2 
72.0 


0.5 
3 
1 
5 
0.7 
3.5 


1000 
1000 
1000 


III 


Corresponds to two standard deviations from 2loud -center concentration. 


bAfter H+! hr. 


Table 2. 
Potential surface exposure and downstream exposure rates. 


Fraction 
Distance 
Yield 
deposited 
downwind 
(kt) 
(%) 
(km) 


1 
100 
10 
100 
1000 


10 
10 
10 
100 
1000 
100 
J 
10 
100 
1000 


Infinite 
v/hole-body 
external gamma dose 
(rem) 


External gamma 
exposure rate at 
time of arrival 
(R/hr) 


2 5,000 
200-1200 
0.3-1.0 


5000 
1000-1500 
1-4 
800 
300-350 
1.5-5 


25,000 
35-100 
<0,01 
7000 
150-200 
0.015-0.05 
1200 
45-55 
0.02-0.06 
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Distance from ground zero • 


1000 


• km 


Fig. 1. Initial gamma exposure rate due to the deposition of the vertical integral at 
various distances from ground zero. The upper curve for each yield repre­ 
sents the case of slow horizontal diffusion; the lower curve represents the 
case of fast diffusion. 
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Distance from ground zero — km 


Fie 2 
Vertical integral (infinite whole-body exposure) due to gross gamrnc; radiation 
as a function of distance from ground zero. The upper curve for each yield 
represents the case of slow horizontal diffusion; the lower curve represents 
the case of fast diffusion. 


- 6 


10-kt debris, and the lower 1% of the 
100-kt debris are available for scav­ 
enging. 
With these limits as a basis. Table 2 
summarizes exposures and exposure 
rates for each of the three yields for 
distances of 10, 100, and 1000 km down­ 
stream from ground zero. 
Observe that 
while the total exposures due to the 10- 
and 100-kt clouds are reduced by factors 
of 10 and 100 respectively, rather high 
values can still occur even out to 100 km. 


In the preceding section of this report 
it has been assumed that nearly all the 
debris in a vertical column from the 
ground up to the top of a rain cloud is 
scavenged hv precipitation. 
But the 
presence of rain clouds implies the 
existence of dynamic effects which can 
greatly al'er the vertical distribution of 
debris from that in the cloud-free environ­ 
ment. 


Consider, for example, a simplified 
picture of the dynamics of a growing 
cumulus cloud. 
Moist unsaturated 
environmental air is drawn up through the 
base of the cloud and is accelerated up­ 
ward because of buoyancy caused by the 
release of latent heat of vaporization. 
This indrawn air either remains a part 
of the young growing cloud or it leaves 
through the anvil top of a large fully 
developed cloud. The debris tends to 
move with the air so that the concentration 
of debris at all heights in a cumulus cloud 
developing in this manner would tend to 
be the same as the debris concentration 
of the environment in the vicinity of the 
base of the rain cloud. Now if the debris 


It should be recalled, however, that the 
10% and :% availability of the 10- and 
100-kt cloud debris was based upon a 
specific case and may be greater or less 
depending upon individual situations. 
For 
2 
example, Conner and Phillips 
in a 
similar calculation founa that 100% of 


the radioactivity from the 10-kt cloud 


was available for scavenging based upon 
his calculations of stabilization neipht, 
cloud thickness, and height of scavenging 
systems. 


concentration is maximum near the rain 
cloud base, the amount of debri.^ in a 
vertical column in the rain cloud will be 
much larger than the amount of debris in 
a corresponding column of the clear 
atmosphere. 
When rain falls, mo±e than 
one times the vertical intergral of the 
debris concentration in the clear atmos­ 
phere can be brought to the ground 
beneath the debris-cloud center. 
In addi­ 
tion, the mass of air that passes through 
the rain cloud may be much larger than 
the mass of air in the cloud itself, further 
increasing the amount of debris that can 
be collected by drops and droplets and 
eventually brought to the ground. 
On the 
other hand, if the debris is primarily 
above the base of the rain cloud, it can 
be imagined that the cumulus cloud would 
grow up through the contaminated regions, 
pushing the debris aside and thereby 
preventing collection of debris by rain­ 
drops. 
From these simple sxamples it 
can be seen that the dynamics of individual 
convective cells are a major factor in con­ 
trolling the amount of debris that enters 
the rain cloud and interacts with drops. 


Dynamic Effects of ConvecMcn 
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As debris moves through a rain cloud, 
some of the radioactivity is collected by 
raindrops and cloud droplets. 
Some of 
the cloud droplets are collected by falling 
raindrops v ule others evaporate when 
the cloud dissipates. 
The contaminated 
raindrops eventually fall to the ground, 
carrying with them all of the radioactive 
debris t'.iey have collected. 
Therefore, 
the debris in a rain cloud that is collected 
by raindrops is scavenged, and the radio­ 
activity that does not get into raindrops is 
left in the atmosphere. 
Consequently the 
microphysical processes by which debris 
particles get into raindrops are the major 
factors that determine what fraction of 
the debris passing through a precipitating - 
system is actually scavenged. 


The remainder of this section is devoted 
to a description of a simple model of 
cumulus convection that has been used to 
estimate the typical range of values for 
the ratio of the vertical integral of the 
debris concentration inside a rain cloud 
to a similar vertical integral in the cloud- 
free environment. 
This ratio of vertical 
integrals is given the name enhancement 
factor, although enhancement factors 
less than 1 correspond to cases in which 
the rain cloud contains less debris than 
the environment. 
In following sections 
of this report, the scavengabiiity of debris 
iarticlos and the microphysical processes 
by which debris is collected and brought 
to the ground with falling raindrops will 
be discussed. 


DESCRIPTION OF THE MODEL 


A one-dimensional steady-state cumulus 
dynamics model has been used to estimate 
the importance of the dynamics of cumulus 
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convection on the scavenging of dsbris 
from the atmosphere. 
This model was 
developed by Weinstein and Davis (1968). 
It follows a sample parcel of air rising 
through a rain cloud; the cloud as a whole 
is assuned to be made up of many of 
these parcels one on top of the other and 
each with the properties of the sample 
parcel at the appropriate level; hence the 
structure of the rain cloud does not change 
with time. 
The sample parcel is initially 
at cloud base and has a small upward 
velocity. 
The surrounding environment 
has a vertical profile of temperature and 
moisture content determined from a 
standard radiosonde sounding. 
Calcula­ 
tions of the cloud tempe ature, vertical 
velocity, mixing ratio, cloud water content, 
hydrometeor water content, and undraft 
radius are mado for the sample parcel at 
each level in the model cumulus cloud. 
The sample parcel also interacts with the 
surrounding environment by entrainmem, 
The parameterized entrainment rate, \x. 
is taken from experimental studies of 
buoyant jets and bubbles and is given by 
the expression 


1 dM 
0.2 


M 
M dz 
R " 
c 


where M is the mass of air, z is height, 
and R 
is the radius of the core of the 
updraft. 
Further details of the dynamic 
and microphysical model are given in the 
above cited reference. 
fn addition to the rain cloud properties 
already mentioned, calculations of the 
concentration per unit mass of debris at 
each level in the cumulus cloud are in­ 
cluded in the numerical model. 
The only 
process that alters the concentration of 


debris within the sample parcel is e n t r a p ­ 
ment of environmental air with a differ­ 
ent debris concentration; therefore th? 
concentration per unit mass of debris in 
the cloud at level t, x^, can be related 
to the concentration per unit mass of 
debris in the cloud at the previous lower 
level, Xn,j' hy the expression. 


x l ' 
1 + /J Az 
' 


where C„ is the concentration per unit 
mass of debris in the environment outside 
the cloud at level i, y is the mass e n t r a p ­ 
ment rate, and Az is the thickness of the 
layer. 
At the lowest level of the model 
(at base of the cumulus cloud) the debris 
concentration is assumed to be equal to 
the environmental debris concentration. 
The numerical value of the entrainment 
rate, M, depends on the updraft radius 
of the cloud which is, in turn, dependent 
upon the other physical properties of the 
rain cloud and the environmental sounding. 
Therefore, each cumulus cloud has a 
unique vertical distribution of debris just 
as each r3in cloud has its own vertical 
<JistviWt\<in of drop concentrations, vrp- 
draft radius, vertical velocity, etc. 
The 
question of whether debris is collected by 
drops or droplets is not considered in 
this simple model although it is of crucial 
importance when considering the scaveng­ 
ing of debris. 


DESCRIPTION OF THE PRECIPITATING 
CLOUD 


The radiosonde sounding for Salt Lake 
City on 16 May 1971 at 1200Z was used to 
describe the environment in which pre­ 
cipitating clouds developed. 
Scattered 


showers were occurring in the area at 
this time. 
The radar-measured tops of 
rain clouds were between 6000 m and 
9000 m, and their bases were at 2700 m. 
The assumed updraft radius at cloud base 
was 1 km and the assumed updraft velocity 
at cloud base was 2 m sec 
. 
The 
freezing temperature for ail drops and 
droplets was assumed to be -25°C, a 
temperature which occurred at an 
altitude of about 6100 m. 
The thick­ 
ness of the layers in the model was 
200 m. 


The environmental and in-cloud tem­ 
perature profiles are displayed in Fig. 3. 
At low levels the in-cloud temperature 
is higher than the environmental tempera­ 
ture, 
primarily because of the release 
of latent heat as water vapor condenses. 
This positive buoyancy accelerates the 
parcel upwards as can be seen in the plot 
of vertical velocity in Fig. 4. 
The jog 
in the in-cloud temperature profile at 
6100 m is due to the sudden freezing of 
all drops and droplets at this level and 
the resultant release of the latent heat of 
fusion. 
This positive buoyancy leads to 
a corresponding increase in vertical 
velocity at levels j-ast above 6100 m and 
indicates the importance of freezing in 
the dynamics of cumulus convection. 
Above 7300 m the buoyancy of the parcel 
is negative and the vertical velocity drops 
to zero at 7900 m, which is, by definition, 
the top of the rain cloud. 
Two other 
factors also tend to decrease vertical 
velocity: the drag of drops and droplets 
which are carried along with the parcel, 
and entrainment of stationary environ­ 
mental air. 
A cloud top of 7900 m is 
consistent with radar measurements of 
the heights of rain clouds in the Salt Lake 
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Fig. 3. Variation of temperature with height in the environment and in the model 
cloud. 
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Fig. 4. 
Distribution of vertical velocity with height In the model cloud. 


City area at the time of the radiosonde 
sounding. 
The vertical distributions of cloud 
droplet mixing ratio and raindrop mixing 
ratio of hydrometeor water are shown in 
Pig. 5. Cloud droplets are formed by con­ 
densation of water vapor in the rain cloud 
and are removed by collection by hydro- 
meteors. Raindrops are formed from 
cloud droplets and grow by coalescence. 
There is no mechanism in the model for 
a decrease of hydrometeor water in the 
sample parcel; therefore the raindrop 
water content must increase with height 
and be a maximum at the top of the cloud, 
in natural rain, fallout of drops would 
greatly alter this particular profile. 
The 
discontinuity in the slope of the cloud 
droplet concentration curve at 6100 m is 
due to the freezing of drops at that level, 
and to changes in certain microphys- 
ical parameters at the freezing level. 
If all the hydrometeor water in the 
rain cloud at any instant fell to the 
ground, there would be 1.6 mm of 
rainfall. 


The characteristics of cumulus clouds 
that affect the vertical distribution of 
debris most directly are the updraft 
radius and the entrainment rate; these 
quantities are plotted in Fig. 6. 
The 
mixing rate is greatest in the region of 
smallest updraft radius which corresponds 
to the height of the largest updrafts. 
At 
levels where the vertical velocity is 
small the entrainment rate is also small. 
At the top of the rain cloud approx­ 
imately two-thirds of the mass in the 
sample parcel is air that has been 
entrained by the rising parcel; the 
other one-third is air that entered at 
cloud base. 


- 1 1 - 


DEBRIS DISTRIBUTIONS 


The environmental debris distribution 
that was used for the model calculations 
was taken from results of a 2BPUFF 
diffusion code (Crawford, 1966) 
calcula­ 
tion. 
The debris concentration 100 km 
downwind from the detonation in the free 
atmosphere of a 1-kt all-fission device 
has a Gaussian profile with a peak con­ 
centration of 5 X 10~ 
3 Ci/m 
3. 
The height 
of this peak concentration is 2800 m and 
the standard deviation of the profile in 
the vertical is 500 m. 
It is assumed that 
the variations of concentration in the 
horizontal can be neglected over the 
region that interacts with the model 
cumulus cloud, and values at the center 
of the debris distribution in the horizontal 
are used. 


The vertical profiles of debris con­ 
centration in the environment and in the 
rain cloud for a 1-kt free-air burst which 
interacts with the model cloud after 
moving downwind 100 km are shown in 
Fig. 7. 
Since die peak environmental 
debris concentration occurs at the base 
of the rain cloud, the debris concentration 
in the cloud is large. 
Mixing of cleaner 
environmental air as the parcel rises 
c 
reases the debris concentration in the 
cloud, but because entrainment is a 
relatively inefficient process the debris 
concentration remains high at all levels. 
Since the cumulus cloud model does 
nothing to alter the debris concentration 
below the rain cloud base, both profiles 
are identical below 2700 m. 
The enhance­ 
ment factor, which has been defined as 
the ratio of the vertical integrals of the 
two profiles in Fig. 7, has a value of 
3.04. 
This implies that if all the debris 
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Fig. 5. Vertical distribution of liquid water components (cloud droplets and rain drops) 
in the model cloud. 
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Fig. 6. Variation of entrainment factor and size of the updraft with height for the 
model cloud. 


-12- 


3 
3 
Debris concentration — (xlO 
Ci/m ) 
Fig. 7. Vertical distribution of debris in the environment and in the model rloud for 
a 1-kt all-fission free-Ear burst. 
The debris is assumed to interact with the 
rain cloud after moving downwind 100 km with the rain cloud where the height 
of the peak concentration for the Gaussian environmental profile is at 2.8 km 
and the standard deviation in the vertical is 0.5 km. 


in and below this cumulus cloud were 
brought scavenged, approximately three 
times the vertical integral of debris in 
the environment outside the c, mulus 
cloud would be brought to the ground with 
the rain. 
It is apparent from the previous example 
that the d oris concentration at the base 
of the r- m cloud is a very important factor 
in controlling the debris concentration in 
the entire rain cloud. 
Therefore, addi­ 
tional sample calculations hr"° been 
performed for several other environ­ 
mental debris distributions interacting 
with the same cumulus cloud; the only 
change is in the height of the peak debris 


concentration. 
These results are 
illustrated in Figs. 8 through 12 and are 
summarized by the data in Table 3. 
From the data it can be seen that the 
enhancement factor is largest when the 
maximum value of concentration occurs 
at the base of the rain cloud. 
When the 
peak concentration is below cloud base 
the enhancement factor is still greater 
than 1, except for the limiting case when 
all the debris is below cloud base, then 
the two profiles are identical and the 
enhancement factor is 1.0. 


When the height of che peak environ­ 
mental concentration is found at greater 
levels above the base of the rain cloud. 
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Debris concentration 
— radioactivity volume 


Fig. 8, 
Vertical distribution of debris in the environment and in the model cloud for 
a hypothetical free-air burst with the assumed height of the peak concentra­ 
tion of the Gaussian environmental profile at 2.2 km and the standard devia­ 
tion in the vertical of 0.5 kin. Numerical values for the ordinate would depend 
on the total amount of radioactivity released by the free-air burst. 


Debris concentration — radioactivity/volume 


Fig. 9. Same as Fig. 8, except that the assumed height of the peak environmeiital 
concentration is 1.6 km. 
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Debris concentration — radioactivity/volume 


Fig. 10. 
Same as Fig. 8, except that the assamed height of the peak environmental 
concentration is 3.4 km. 


1 
2 
3 
4 


Debris concentration — radioactivity/volume 


Fig. 11, Same as Fig. 8, except that the assumed height of the peak environmental 
concentration is 4.4 km. 


2 
3 


Debris concentration — radioactivity/volume 


Fig. 12. 
Same as Fig. 8, except that the assumed height of the peak environmental 
concentration is 5.0 km. 


Table 3. Important parameters for comparison of debris profiles. 


Figure 


Height of peak 
environmental 
concentration 
(m) 


Standard deviation 
in the vertical of 
environmental 
concentration 
(m) 


Concentration 
at cloud base 


(X 10" 
3 Ci m~ 
3) 
Enhancement 
factor 


7 
2800 
500 
4.90 
3.04 


8 
2200 
500 
3.03 
2.32 


9 
1600 
500 
0.44 
1.20 


10 
3400 
500 
1.88 
1.60 


11 
4400 
500 
0.02 
0.73 


12 
5000 
500 
0 
0.60 


13 
4100 
1000 
1.88 
1.21 


14 
4500 
1000 
0.99 
0.88 


transport of debris into the cumulus cell 
by entrainment becomes more important, 
since only a small amount of debris 


enters through cloud base. 
At levels 
where the debris concentration in the 
environment is greater than the debris 
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concentration in the rain cloud, e n t r a p ­ 
ment increases the debris concentration 
in the cloud. 
Conversely, at levels where 
the debris concentration in the environ­ 
ment is smaller than the debris concen­ 
tration in the rain cloud, entrainment 
decreases the concentration of debris in 
the cloud. 
Consequently, the peak con­ 
centration in the precipitating cell occurs 
above its base and also above the height 
of the peak environmental concentration. 
Also, the concentration of debris in the 
rain cloud in lower regions is smaller 
than the concentration of debris in the 
environment, but at higher levels, once 
again, the precipitation region contains 
more debris. 
The value of the enhance­ 
ment factor, which is maximum when the 


environmental debris concentration has 
its peak value at the base of the rain cloud, 
decreases as the height of the peak 
environmental concentration increases. 
The minimum value of the enhancement 
factor of 0 when the debris is completely 
above the rain cloud. 
The enhancement 
factor is greater than 1 only when signif­ 
icant amounts of debris enter the pre­ 
cipitating system through its base. 


Additional representative calculations 
of debris profiles for the same model 
cumulus cloud are illustrated in Figs. 13 
and 14; here the standard deviation of the 
environmental profile has been doubled. 
These results can be compared with two 
cases of the previous results, one for 
which the concentrations at the base of the 
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Debris concentration — radioactivity/volume 


Fig. 13. 
Same as Fig. 8, except that the assumed height of the peak environmental 
profile is 4.5 km and the Gaussian standard deviation is 1.0 km. 
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2 
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Debris concentration — radioactivity/volume 


Fig, 14. 
Same as Fig, 8, except that the assumed height of the peak environmental 
profile is 4.1 km and the Gaussian standard deviation is 1.0 km. 


rain cloud are the same (Figs. 10 and 13) 
in which case the enhancement factor is 
smaller for the broader profile, and the 
other for which the peak concentrations 
are at the same height (Figs. 11 and 14). 
in which case the enhancement factor is 
larger for the broader profile. 
This 
latter result follows primarily from the 
fact that the value of concentration at the 
base of the precipitation system is much 
larger for the broader environmental 
profile. 


In summary, it has been shown that 
values for the enhancement factor will 


typically range from 0 to values of 3 or 
4. 
The highest values occur when the base 
of a rain cloud forms near the level of 
peak concentration, in which case the 
amount of debris available for scavenging 
can be several times ss large as the 
amount of debris that would have been 
present if the dynamics of convention had 
been ignored. 
When the debris is mostly 
above the base of the rain clouds the 
enhancement factor is less than 1 imply­ 
ing that the dynamics of convection tend 
to reduce the amount of debris that can 
be scavenged. 
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A Microphysical Description of Precipitation Scavenging 


Aerosol particles are removed from 
the atmosphere by different processes 
depending on their size. 
Particles 
larger than about 1 ^m can be scavenged 
by direct impaction with raindrops. 
This 
inertial capture process is the primary 
removal mechanism for washout or below- 
cloud scavenging, the washout efficiency 
being a function of particle size. 


For rainout or in-cloud scavenging, 
other processes must be considered in 
addition to inertial capture. 
Aerosol 
particles with perfectly wettable surfaces 


and radii greater than 0.1 ium can serve 
as nucleation sites for the formation of 
cloud droplets. 
Soluble particles with 


radii as small as 0.01 ^m can be 
nucleated in clouds. 
According to argu- 
ments by Fletcher 
and experimental r e - 
4 
suits by Twomey , an insoluble particle 
whose surface departs from being com­ 
pletely wettable, i.e., the contact angle 
is greater than say 6°, cannot serve as a 
condensation nucleus for typical cloud 
supersaturations. 
Nucleated particles 
can be deposited on the ground either 
through growth by diffusion to raindrop 
size and falling (a very rare event) or 
growing to several microns in radius and 
being accreted by a raindrop. 


Particles of radius less than about 
0.02 /jm display considerable Brownian 
motion under atmospheric conditions, 
allowing them to collide and attach to drop­ 
lets or raindrops. 
This scavenging 
mechanism is directly effective only for 
short times after the particles are formed 
since, given time, they will readily attach 
to larger debris particles or to natural 
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aerosol particles. 
Those po.rt 
;cles 
attached to droplets due to Brownian 
capture can be 2-emoved to the ground 
through accretion of the droplet. 


Fallout particles from free air bursts 
consist mainly of oxides of the dominant 
casing materials. 
The various radio­ 
nuclides make up a small percent of the 
total particle mass. 
Indications are tt- it 
particles are insoluble in water and are 
wettable due to their large surface energy 
compared to that of water. 
Some of the 
radionuclides are sc'-ible although it is 
not clear if the amount of soluble material 
on the particle surface can significantly 
change the nucleation characteristics as 
suggested by Hicks. 


Indications are that the particle size 
resulting from free air bursts are distrib­ 
uted log-normally and that the mean particle 
13 14 
size decreases with increasing yield. 
' 


It appears reasonable to expect low yield 
air bursts to produce particles suscep­ 
tible to inertial capture and nucleation 
scavenging. 
However, higher yields may 
produce particles below the nucleation 
threshold, leaving Brownian scavenging as 
the dominant removal mechanism. 


In the remainder of this section an 
attempt is made to formulate th» dominant 
scavenging processes discussed above. 
The model will be concerned with 
describing the removal of debris from a 
parcel of air in which the rain rate is 
constant throughout. 
It is further assumed 
that the concentration of debris particles, 
cloud droplets, and raindrops can be con­ 
sidered spatitlly uniform inside this 
parcel. 


Each of the dominant removal mechanisms 
involve the collision and subsequent col­ 
lection of a debris particle, or a water 
droplet which contains the particle, by a 
raindrop. 
Denoting as n(r, t) dr the 
number of these target spheres of radius 
between r and r + dr per unit volume of air 
at time t, the rate at which such spheres 
are removed by falling raindrops is: 


i ^ f i i = -A{r,t)n(r,t) 
(1) 


Mr, t) is the late at which one target 
particle of radius r is removed from the 
volume at time t. 
A solution to this 
equation can formally be written as: 


n(r, t> = n(r, o) exp | -J 
dt'Mr, t')\, 
(2) 


where 


n(r, o) is n(r, t) evaluated at t = 0. 


\n expression for Mr, t) as given by 
15 
Chamberlain 
can be written as: 


Mr.t) A 
Jo 
dRN(R.t) jrlT[V(R) - v(r)I 


X E(R, r, f). 
(3) 


Here N(R, t) is the number of raindrops 
per unit volume with radius R and fall 
velocity V(R). 
E(R, r, f) is the collision 
efficiency of the raindrop with a target 
sphere of radius r and mass density ?. 
v(r) is the fall velocity of the target sphere. 


Scavenging by diffusiophoresis and 
thermophoresis does not appear to be 
important. 
The thermophoretic force on 
a submicronic particle juxtaposed with a 


droplet u' der'oing diffusional growth 
16 
dominates the diffusiophoretic force 
causing the particle to be repelled from 
the droplet. 
Although the attractive dif- 


17 


fusionhoretic force dominates for particle 
sizes greater than about 1 urn, the r e ­ 
sulting increase in the collision efficiency 
when the diffusiophoretic force is 
accounted for is of little significance. 
The removal rate given by Eq. (3) de­ 
pends on time through the time dependence 
of the drop spectra. 
This expresses the 
fact that throughout the history of a cloud 
the number of drops of radius R in any 
unit volume changes with time, 
Under the 
assumption that the drop size distribution 
is time independent, Eq. (2) reduces to: 


n(r, t) = n(r, o) exp (-A(r)tJ. 
(4) 


For computation of the removal rate, 
Mr), 


the terminal velocities of the drop given 
1 a 
by LeClair, et al. 
are used for Reynolds 
number less than 400 and for larger drops 
19 
those of Gunn and Kinzer. 
Terminal 
velocities of the collected sphere are 
assumed to obey Stoke's law and the col­ 
lision efficiencies used a r e those of 
20 
Langmuir. 
The drop spectra is taken 
21 
from results by Best 
and can be written 


N{R) 
12n W <2R)" 


exp um^y 
(5) 


Here n = 2.25, W is ths mixing ratio and 
a is a constant dependent upon rain rate. 
Values of W and a quoted by Best * for 
various rain rates, I, are given in 
Table 4. 
The drop spectrum is defined 
here such that N(R) dR represents the 
number of drops per unit volume with 
radius between R and R + dR. 


Debris particles of radius greater 
than about 1 urn are subject tc inert ial 


-20- 


Table 4. 
Values of a and W correspond­ 
ing to different rain rates, I. 
From Ref. 21. 


T 
mm 
1 ~ET 
a(cm) 
WX 10 
9 (cm/cm) 
3 


0.5 
0.1107 
37 


1.0 
0.13 
67 


2.5 
0.1608 
145 


5.0 
0.1888 
261 


10.0 
0.2218 
470 


25.0 
0.2743 
1020 


capture by raindrops. Mathematically this 
situation can be described by Eqs. U)-t4) 
where the collected particle is taken as a 
spherical debris particle. 
From Eq. (4) 
the fraction of debris particles of radius r 
which have been removed from the parcel 
at time t is: 


1 ' 
n(r, o) 


1 - expf-Aj 
Mil. 
(6) 


Using values of the drop spectra found by 
22 
Mason and Andrews 
for warm frontal 
showers, the washout coefficients are about 
a factor of 2 smaller than those shown in 
Fig. 15. This indicates that the washout 
coefficient depends on type of rain shower 


as well as the rain rate. 


In the present model, nucleation 
scavenging is calculated by assuming every 
debris particle of radius 0.1 m <_ r < 10 ^m 
serves as a condensation nucleus. 
The 
resulting droplet is assumed to immediately 
grow to 10 jum radius. 
Since the mass 
density of debris particles is typically 
greater than the density of water and of 
most natural aerosol particles, and 
further since there may be a significant 
number of debris particles with a radius 
of several microns, the density of the 
debris nucleated droplet is taken as: 


? = S r 
3 + £ 
P 
* w W"')} "V 
"> 


Assuming that every particle captured by 
a drop gets deposited on the ground, 


Fjfr, t) gives the fraction of particles of 
radius r origfriaYi-y ir the air that react 
the ground as a function of time. 
Considering tnertial capture to be the 
only scavenging mechanism for washout, 


Eq. (6) gives the fraction of particles of 
radius r removed from a below-cloud 
parcel as a function of time. 
Figure 15 
gives calculated values for the washout 


efficient, A^r), as a function of particle 
size for various rain rates assuming a 
3 
particle mass density of b g/cm , 
These 
values substantually agree with those found 


by Chamberlain. 
Values of the washout coefficient are 
sensitive to the s i t e distribution, of drojis. 
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Here ? is the density of the debris parti­ 
cle of radius r and f 
is the density of 
water. 
The radius of the nucleated drop­ 
let is r w . 
The fraction oj debris particles of 
radius r removed from the parcel by 
nucleation and subsequent impaction with 
a raindrop is given in analogy with Eq. (5) 
and under the above assumptions as: 


F 2(r,t) = ! - exp (-A^rlt). 
(8) 


A„(r) is the removal rate for particles of 
radius 10 Mm and density as shown in 
Eq. (7). 
Figure 16 represents the rainout 
coefficient, A,(r), f o r 0 , l u m < 
r < 10 ^m. 


For particles of radius greater than 10 urn, 
the rainout and washout coefficients are 
3 
equal sise since K, - % = 5 g/cm . 
For 


Rain rate = 25 mm/hr 


10 mm/hr 


5 mm/hr 


2.5 mm/hr 


1 mm/hr 


J 
i 
L 
1.2 
1.4 
1.6 


Particle radius — (xlO 
cm) 


Fig. 15. Washout coefficient versus particle radius for various rain rates. 


a given rain rate. Fig. 16 shows the rain- 
out coefficient to be relatively constant 
for the particle sizes considered. 
In order to make a complete rainout 
assessment, one should have details of 
the size distribution of debris particles. 
However, as a result of the assumption 
that every debris particle between 0.1 
and 10 /jm radius serves as a condensa­ 


tion nucleus when subjected to a cloud 
environment, a further simplification 
concerning the rainout calculation becomes 
apparent. Let m(r) dr represent the 
mass contained in debris particles of 
radius between r and r + dr. The total 
mass of debris deposited on the ground by 
particles of radius greater than 0.1 jjm 
becomes: 
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Fig. 16. 
Rainout coefficient versus particle radius for various rain rates. 
1 
M(t) = I 
dr m(r) 
0.x ^m 


M(t» 


X [1 - exp (-^(rjt)!, 
(9) 
(l - e " 
A 2 ) M o. 


dr m(r) 


(10) 


Here M is the mass contained in particles 
o 
r 
Here A„ is the rainout coefficient for a 
particular rain rate as plotted in Fig. 16. 
of radius greater than 0.1 p. Under the 
Since A„(r) is a weak function of r, Eq. (9) 
further assumption that the radioactivity 
can be approximated by: 
is contained in particles of greater than 
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0.1 urn radius, which is approximately 
An evaluation of the potential hazard 
true for yields as low as 1 kt, the fraction 
due to rainout from a 1-kt free-air burst 
of radioactivity deposited on the ground as 
based on the above results can now be 
a function of time is plotted in Fig. 17. 
made. Consider the debris system to 


Time — min 
Fig. 17. Fraction of radioactivity originally in the rain cloud deposited on the ground 
versus time for various rain rates. 
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Fig. 18. 


o 25 mm/hr rain rate 


a 2.5 mm/hr roin rate 


1000 
Distance from ground zero — km 
Data points represent initial gamma exposure rates due to rainout of debris 
from a 1-kt free-air burst versus distance from ground zero. Rain is 
assumed to fall for 15 min with rain rates of 2.5 and 25 mm^r. The solid 
curves are reproduced from Fig. 1. 
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move downwind and to diffuse according 
to the Lagrangian diffusion code (2BPUFF) 
under non-shear conditions with horizontal 
dissipation, e = 1.8 ergs/g-sec. At vari­ 
ous distances downwind the debris system 
is assumed to interact with a cumulus 
cloud, as described previously, giving an 
enhancement factor of 3. The radioactive 
particles are assumed uniformally 
distributed within the rain cloud and 
deposited by rainout as a function of time 
and for various rain rates as given in 
Fig. 17. Washout is neglected and it is 
assumed that the drops fall on the area 
directly below the rain cloud. The data 
points in Fig. 18 show values of the initial 
infinite whole body exposure rate versus 
distance downwind for rain rates of 25 and 
2.5 mm/hr lasting for 15 min. The curves 
in Fig. 18 are reproduced from Fig. 1. 
These new results on the potential 
hazards from rainout of debris from the 
1-kt cloud are worthy of special comment, 
namely, that the results of the present 
analysis indicate that the potential dose as 
a function of range for rainout (Fig. 18) 
constitutes a reasonable estimate of the 
expected hazard for the conditions cited. 
This expected estimate level of hazard 
is quite comparable to that previously 
reported in "Potential Exposures from 


Low-Yield Free Air Bursts." 
1 During 
the remainder of this contract year, the 
LLL Rainout Team will continue to refine 
the present estimate, as represented in 
Fig. 18, for effects of wind shear (not 
expected to be large in summer), and 
other assumptions inherent in the analy­ 
sis. Further work is required in regard 
to the variation of particle size distribu­ 
tion with yield in order to extend the 
present analysis techniques to yields other 
than 1 kt. This work will be a portion of 
our proposal for the coming year. In the 
future we will also evaluate the probability 
associated with meteorological conditions 
consistent with Fig. 18 and other scenarios 
of interest to DNA. 
Rainout estimates assuming nucleation 
scavenging and inertial capture were 
possible for a 1-kt free-air burst be­ 
cause in this case it appears that very few 
o" the debris particles are smaller than 
23 
0.1 ftm radius. 
It appears that the mean 
particle size decreases with increased 
yield such that for higher yields a 
significant portion of the radioactivity is 
contained in particles smaller than0.1 jum. 
This being the case, different scavenging 
mechanisms (probably Brownian capture) 
become important and the above analysis 
does not apply. 


Thoughts on Altering the Scavenging of Debris Particles 


The problem here is to change or to 
control the rate of rainout from clouds. 
We are studying the possibilities of 
changing their macrophysical dynamics 
by interfering with their microphysical 
processes. One can add to a cloud or 
take parts of it away, and at first look 
it appears that adding something is the 
easier operation. 


There are a number of possible effects 
of additives on the microphysical proc­ 
esses of water and ice particle formation 
in a cloud that can influence the overall 
macrophysical dynamics of the cloud. 
Additives in the form of an excess num­ 
ber of nucleating agents could hinder the 
development of water and/or ice particles 
to a size large enough to precipitate out 


26- 


of a cloud. 
Small amounts, a monolayer 
or even a partial monolayer, of materials 
with low surface energies adsorbed on to 
the high energy surfaces of nucleating 
particles can change their wetting charac­ 
teristics and thus their nucleating rates 
and possibly put them into the nonnucleat- 
ing category as far as the cloud is con­ 
cerned. 
Small amounts of materials can 
also affect the electrical properties of the 
water and/or ice particles, the elastic 
and inelastic collision processes, the 
accommodation and reflection coefficients 
for the water molecule on water and ice, 
and change the condensation and/or evap­ 
oration rates of water molecules between 
water droplets and between water and ice 
particles. 


Let us look at the use of additives to 
modify the rate of nucleation of particles 
by changing their wetting characteristics. 
We shall assume the validity of the 
Fletcher analysis 
of heterogeneous nucle­ 
ation of water vapor by insoluble particles 
as a function of particle size, contact 
angle and saturation ratio in the cloud. 
Figure 19 is a set of Fletcher's curves 
for a given nucleation rate. 
Since the 
saturation ratio in most clouds is about 
1.001 and only under unusual circum­ 
stances may reach as high as 1.01, it can 
be seen from Fig. 19 that a particle with 
radius less than 10 
cm will not act as a 


nucleation center. 
Even particles with a 
-5 
radius greater than 10 
cm need to be 
almost completely wettable before they 
can act as nucleation centers for water 
vapor in a typical cloud. 


The solid particles that are expected 
to be picked up or injected into a cloud 
would be metal oxides in the form of crys­ 
tals or glasses and possibly some free 


metallic particles. All such particles 
have a large surface energy compared to 
water (73 ergs/cm ). See Tables 5 and 6. 
One would expect water to wet and spread 
29-31 
on all such particles, 
and from the 
Fletcher analysis all such particles would 
be active nucleating agents in a cloud if 
-5 
3 
they were 10 
cm in radius or larger. 
The rates of nucleation by these solid 
high surface energy particles can be 
changed by the addition of low surface 
energy materials that will adsorb onto 
the surface of these particles in the cloud 
and lower their surface energies to less 
than that of water. 
A monolayer of 


Table 5. Surface energies of various 
materials, taken from refer­ 
ences cited. 


7 at 0°K 
Materials 
(erg/cm 
2) 
Reference 


Oxides of General Formula MO 


MgO 
1090 
24 


F e O 
1060 
24 


MnO 
1010 
24 


CaO 
820 
24 


SrO 
700 
24 


BaO 
604 
24 


BeO 
>1420 
24 


CdO 
530 
24 


LnO 
600 
24 


P b O 
250 
24 


General Formula M „ 0 ? 


B 2 ° 3 
79.5 <§• 900°C 
25 


A 1 2 ° 3 
905 
@ 1850°C 
26 


General Formula MO„ 
uo2 
642 ± 20% 
24 


Z r 0 2 
800 ± 20% 
24 
T h 0 2 
530 ± 20% 
24 
S i 0 2 
>800 
24 
TiQ 2 
>800 
24 
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Pig. 19. Critical supersaturation ratio for nucleation of a water droplet in 1 sec on a 
particle of given radius and surface properties defined bym = cos 0 where $ 
is the contact angle (from Ref. 3). 


perfluoro fatty acid on platinum reduces 
its specific surface free energy from 
2 
about 2000 ergs/cm to about 10 ergs/ 
2 32 33 
cm . 
' 
The water should not wet nor 
spread on these particles. Even a small 
departure from perfect wettability will 
hinder the nucleation dynamics. There 
are a number of possible materials that 
could be selected, and these could be 
introduced into the cloud by rockets, bal­ 
loons, smudge fires, antiaircraft guns, 
the device, etc. If one could make the 
additive available to be adsorbed onto the 
particles not much material would be 


needed. From Table 7 it can be seen that 
3 
for each cm of metal or metal oxide, 
subdivided into 10 
particles with sides 
-5 
-5 
of 10 
cm, it would take only 10 
mole 
of adsorbable material to give all parti­ 
cles a monomolecular coat. 
For optimum control of rainout one 
may want to use one or more types of 
additives at the same time. It may be 
best to keep the insoluble particles from 
wetting by an adsorbable additive and to 
overnucleate the cloud with NaCl parti­ 
cles with masses equal to or greater than 
_2 
10 
grams per particle. This would 
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keep the insoluble particles from nucleat­ 
ing and furnish enough water soluble par- 


Table 6. Surface energy of some metal­ 
lic elements (at their melting 
p o i n t s ) . 
2 7 - 
2 8 


Metal 


Melting 
point 
CO 


Surface- 
energy, 7 
(ergs/cm 
2) 


Aluminum 
660 
900 


Antimony 
630 
370 


Beryllium 
1400 
1000 


Cadmium 
321 
620 


Cobalt 
1495 
1530 
Copper 
1083 
1100 
Gold 
1063 
1120 
Iron 
1534 
1500 


Lead 
325 
450 
Magnesium 
650 
560 
Mercury 
-39 
460 
Nickel 
1453 
1700 


Platinum 
1769 
1800 
Silver 
961 
920 
Tin 
232 
570 
Tungsten 
3410 
2300 
Uranium 
Zinc 
420 
790 


Table 7. Surface area of a 


Length of side 
Number of 
(cm) 
cubes 


10° 
10° 


IO- 
1 
10 
3 


i o - 
2 
10 
6 


10-3 
10 
9 


i c f 
4 i [= 1 (/m) 
io 
1 ; 


IC" 
5 
10 
1 1 


io- 
6 
10 
1 1 


io- 
7 
10 
2 


to" 
8. h 1 A) 


We 0assume that it takes 10 
molecules 
6 X lf>"'' molecules. 


t i d e s for nucleating centers to keep the 
droplets of water and particles of ice 
small enough so that they will not precipi­ 
tate out and carry out the solid particles 
by impact. 


There are a number of possible effects 
of additives on cloud dynamics about which 
little or nothing is known. Research into 
these effects is desperately needed. 
It should also be kept in mind that a 
cloud is a dynamic system, and the appli­ 
cability of equilibrium arguments is some­ 
what limited and may at times be mislead­ 
ing. The Knudsen equation, for instance, 
gives too large a gas transport to and 
from a water drop under nonequilibrium 
conditions. It is probable that the con­ 
densation and evaporation coefficients of 
water molecules at liquid water surfaces 
can be modified by as much as a factor 
of 10 by contaminating the water sur- 
faces 
or even aging them. 
There 


are many different combinations of fac­ 
tors involved in the interaction between 
incident water molecules and surfaces of 
water drops and ice particles. 
These 
factors are the: 


1-cm cube with subdivision. 


Total surface 
Moles of material 
area (cm 
2) 
to cover surface 
3 . 


6 X 10° 
10" 
1 0 


6X 10 
1 
io- 
9 


6 X 10 
2 
lO' 
8 


6 X 10 
3 
io- 
7 


6X 10 
4 
lO' 
6 


6 X 10 
5 
lO" 
5 


6X 10 
6 
lO" 
4 


tSX 10 
7 
lO- 
3 


cover each em 
and each mole has 
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• 
specific condensation coefficient, 
• 
specific evaporation coefficient, 
• 
thermal accommodation coefficient, 
• 
momentum accommodation coef­ 
ficient, 
• 
extent of compliance with the cosine 
law for diffuse reflection, and 


• 
extent the cosine law would obtain 
due to scattering from surface 
irregularities alone. 
37 
Even though Young 
gave a correct 
form of the interaction of a liquid at equi­ 
librium on a solid surface in 1805 it was 
38 
only in 1959 that Johnson 
gave a clear 
statement of the problem and a thermo­ 
dynamic justification of Young's equation 
along the lines of Gibbs' derivation. The 
notion of a wetting tension was given by 
Dupre (1869). 
3 9 and Cooper and Nuttall 
(1915) 
gave the conditions for spreading 


and nonspreading. 
Harkins (1920-1941) 
further developed the notion M a spread­ 
ing tension more from the point of view of 
a gradient of the free energy of the sys­ 
tem with respect to the area. All of the 
above developments are for equilibrium 
conditions and assume no electrification 
of the system. 
The nWecular and/or 
atomic interaction mechanisms of wetting 
and wetting rates are not well understood. 
If one considers small drops on surfaces, 
ignoring inertial effects, one might expect 
the wetting rate of a surface to be ex­ 
pressed by equating the rate of change of 
the reversible work due to "capillary" 
forces to the rate of energy dissipation 
due to the viscous flow of the film over 
the surface. There are a number of efforts 
needed to understand the off-equilibrium 
dynamics in the microphysics of clouds. 


Dry Deposition 


One method by which radioactivity 
from air burst debris clouds may be de­ 
posited upon the ground, generally over 
long periods of time, is through the process 
of dry deposition. In this process the very 
small particles in the clouds are dispersed 
downward through the atmosphere by dif­ 
fusive eddies. As the small particles 
reach the ground, they become attached to 
various surfaces which they may encoun­ 
ter and, unless removed, become sources 
of radioactivity at these points. 
The 
amount deposited upon a given surface 
depends upon a so-called velocity of dep­ 
osition (or deposition velocity) defined 


as 
2 
Amount deposited per cm 
V = 
of surface per sec 
d 
Volumetric concentration 
per cnv* above surface 


,15 


where the volumetric concentration is 
measured very near the surface. 
Values 
of V. vary from 0.1 to 10 cm/sec with 
most reported values near 1 cm/sec (see 
Ref. 15). Thus, to determine the total 
deposition, D, one needs a time history of 
the volumetric concentration x near the 
surface in order to evaluate the integral 


1 


Of course decay must be accounted for 
when considering radioactive materials. 
When the debris cloud from a free-air 
burst becomes stabilized and moves down­ 
stream, the concentration of radioactivity 
near the ground will usually be small 
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(depending upon height of stabilization) 
until the radioactive particles are dis­ 
persed downward by the action of vertical 
diffusion. 
Since this process can be ex­ 
pected to take some time and the amounts 
deposited are small for low concentra­ 
tions, the potential threat from one or 
even a few bursts is seen to be quite 
small in most cases. 
An exception would 
be when the surface rises rapidly toward 
the higher concentrations near cloud cen­ 
ter as in a ^ange of mountains or high 
hills. 
If, however, a large number of 
low-yield bursts were detonated within a 
limited area, the resulting cloud would be 
large in horizontal extent but would be 
stabilized at a low altitude. 
The radio­ 
active particles would soon become dis­ 
persed to the surface, and, due to the 
large horizontal size of the cloud, high 
surface concentrations could be present 
at a given point for long periods depend­ 
ing upon the time of pt-ssage of the cloud. 
Also, as the cloud aged and became very 
large horizontally, even a small surface 
concentration, when integrated over a 
long period of time may result in large 
values of dry deposition; this effect would 
be further enhanced at points far down­ 


stream if conditions were such that no or 
very little deposition of radioactivity 
occurred earlier. 


The effect of wind velocity, i.e., the 
mean speed at which the debris cloud is 
carried downstream, will also have a 
bearing upon the total deposition since, 
for a given distance from ground zero, 
the mean speed will determine the amount 
of decay, and to some extent, the amount 
of dispersion of the radioactivity within 
the cloud. 


Table 8 presents some calculations 
which demonstrate the effects of stabil­ 
isation height and wind speed upon the 
infinite dose from gross gamma due to 
dry deposition over flat ground at various 
distances downstream from ground zero. 
It is assumed that the clouds are created 
by the detonation of 10 Mt of low yield 
devices within an area bounded by a 
radius of 80 km. 
The calculations were 
made using the computer code 2BPUFF 
programmed with a deposition velocity of 
1 cm/sec. 
Trials 1 through 3 demon­ 
strate the relatively small exposures 
obtained from debris clouds which stabi­ 
lize at heights of 5 to 6 km. 
The 
higher values in Trial 2 are caused 


Table 8. Effects of stabilization height and wind speed on the infinite gross gamma 
dose due to dry deposition. 


Trial 


Cloud 
center 
height 
(km) 


Wind 
speed 
(m/sec) 


Initial 
thickness 
(km) 
Dry deposition, R 
No. 


Cloud 
center 
height 
(km) 


Wind 
speed 
(m/sec) 


Initial 
thickness 
(km) 
200 km 
2000 km 
20,000 km 


1 
6 
20 
4.4 
neg. 
neg. 
0.003 


2 
5 
20 
4.4 
0.002 
0.8 
0.2 


3 
5 
10 
4.4 
0.01 
0.005 
0.01 


4 
2 
20 
4.4 
40 
2 
0.15 


5 
1 
10 
4.4 
100 
2 
0.15 


6 
3 
10 
3.0 
neg. 
neg. 
10 
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by experimental use of nigher values 
of the verticr' diffusivity. 
In Trial 4, 
the cloud center is lowered to 2 km and 
the thickness 
is such that a reasonably 
high concentration is present at the sur­ 
face at the time of stabilization. 
The circum­ 
stances are sufficient to cause large expo­ 
sures close- in with decreasing but high 
exposures farther downstream. When the 
cloud center is lowered to 1 km (Trial 5), 
exposures are higher at all distances 
considered even though thp windspeed is 
halved and the travel time therefore 
doubled with respect to Trial 4. 
Trial 6 
was made to demonstrate the effect of a 
debris cloud in which diffusion and depo­ 
sition are inhibited until the cloud is 
20,000 km downstream, at which time a 
large vertical diffusivity is allowed to 
bring debris to the surface. 
This situa­ 
tion might represent the cloud's passage 


A Three-Dimensional A 
Particle-In-Cell 
A Three-Dimensional Atmospheric Diffusion 
Particle-In-Cell Code {ADPIC) 


ADPIC is a numerical three-dimensional 
Cartesian particle diffusion code, capable 
of calculating the time dependent distribu­ 
tion of air pollutants under many condi­ 
tions, which include strongly distorted 
advection wind fields and calm conditions 
for which Gaussian models are not suit­ 
able. 
Basically the code solves the three- 
dimensional advection-diffusion equation 
in its conservative form (pseudo-velocity 
technique), for a given mass-consistent 
advection field, by finite difference 


'The thickness of a debris cloud in which 
the debris is normally distributed in the 
vertical, as is the case here, is the dis­ 
tance between the two standard deviation 
concentration points; i.e., thickness = 4o z. 


over a large cold dome of air and/or pass­ 
ing over an area in which a low-level 
temperature inversion is present for a 
long period of time. 
Furthermore, the 
increase at 20,000 km is partly due 
to the fact that the terrain was pro­ 
grammed to slope upward toward, but 
not to reach, cloud center, 
thereby 
simulating the effect of a debris cloud 
moving over a range of high hills or 
low mountains. 


More work needs to be done on the 
effects of dry deposition. 
Deposition 
velocities over different surfaces and 
velocities of different particle sizes need 
to be considered, and the inhibitory 
effects of inversions and the enhancement 
effect of strong convection and overturn­ 
ing, as well ac the extent to which sloping 
land increases or decreases deposition 
should be studied. 


approximations in Cartesian coordinates. 
The method is based on the particle-in- 
42-44 
cell technique 
with the pollutant 
concentration represented statistically by 
imbedded Lagrangian marker particles in 
an Eulerian grid. 
The inclusion of 
sources, sinks and gravity terms is 
planned. 
The capability of the grid to 
automatically expand and travel with a 
pollutant cloud, makes the code particu­ 
larly suitable for the study of single puff 
releases, such as a free-air burst. 


THE METHOD 


The pseudo-velocity method consists 
of the following: Given the nonlinear 
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transport-diffusion equation, 


| £ + U • VQ = V . (KVQ) 
(11) 


where Q is a scalar concentration, K the 
diffusion coefficient and U the (given) 
mass consistent wind advection field, we 
can, under the assumption of incompress- 
ibility, replace the V - VQ term by 
V^ • (QU). Upon combining the two diver­ 
gence terms we can rewrite Eq. (11) in its 
conservative (pseudo-velocity) form. 


$*V.[Q(U-§*Q)]- 3Q 
3t 


+ V • <QU') = 0 
(12) 


where U' = U - K/Q VQ is the pseudo 
transport velocity. 


The grid mesh of the code is repre­ 
sented by an Eulerian grid consisting of 
three-dimensional rectangular cells of 
uniform size. 
The concentrations Q are 
defined at the centers of the cells and the 
velocities U, U' and -K/Q VQ are defined 
at the cell corners. 
The locations of the 
particles, which represent the pollutant 
cloud statistically, are defined by their 
individual Lagrangian coordinates within 
the Eulerian fixed grid. 


A time cycle of the code is divided 
into an Eulerian step and a Lagrangian 
step and proceeds as follows: 


1. Eulerian Step: The concentrations, Q, 
given for each cell at the beginning of 
the cycle, are used to calculate the 
diffusivity velocities U"D = -K/Q VQ 
which are then added to the given wind 
advection velocities 0 to yield a pseudo 
velocity U' for each cell corner. 


2. Lagrangian Step: Each marker parti­ 
cle contained in a given cell is trans­ 


ported for one time step At with a 
velocity U' 
wlv'rh is computed from 
the pseudo velocities V at the corners 
of the cell, and the particle coordi­ 
nates X by a volume weighting scheme. 


X (new) = X (old) + U' At. 
(13) 


3. Finally, a new concentration distribu­ 
tion, Q, is calculated from the new 
particle positions, which ends the 
cycle. 
The chief advantage of using such a 
hybrid Eulerian-Lagrangian scheme is 
that the fictitious diffusion inherent in a 
purely Eulerian scheme is eliminated. 
The truncation errors inherent in the 
finite difference algoritnms remain, of 
course, and must be minimized by the 
choice of the time-step. 


The boundary conditions are imple­ 
mented in a layer of boundary cells that 
surround the grid. 
In the case of a grow­ 
ing puff cloud the Eulerian grid will be 
expanded automatically, once a particle 
has moved into any of the boundary cells. 
Topography and inversion layer bounda­ 
ries are represented by a deposition and 
reflection boundary conditions, respec­ 
tively. 
The use of relative coordinates 
allows the grid to follow a pollutant cloud 
along its trajectory. 
The computer time 
required for a 15 X 15 X 15 = 3375 cell 
Eulerian grid with 4000 particles is 
approximately one second per time cycle 
on a CDC 7600.* 


Reference tc a company or product 
name does not imply approval or recom­ 
mendation of the product by the University 
of California or the U. S. Atomic Energy 
Commission to the exclusion of others 
that may be suitable. 
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Future Work 


Future plans include continued re­ 
search into the areas outlined in this 
report. 
To account for shear effects on 
the debris system, a three-dimensional 
particle in cell diffusion code capable of 
calculating the time dependent distribu­ 
tion of debris under distorted advection 
"find fields is being developed and 
validated. 


The dynamics of interaction between a 
debris system and a precipitating system 
determine the amount of debris available 
for scavenging and the time history of the 
wet deposition on the earth's surface. 
Future plans include adapting existing 
one- anr". two-dimensional time-dependent 
rain cloud models to describe the interac­ 
tion of a precipitating system with a 
debris cloud. Incorporation of a micro- 
physical description of scavenging into 
these numerical rain cloud models will 
enable estimates to be made of the fate of 
debris within the cloud as a function of 
time and of the spatial and temporal dis­ 
tribution of radioactivity lieposited on the 
ground. 


Theoretical investigations into the for­ 
mation of debris particles to determine 
the particle size distribution for various 
yields and its evolution with time are 
planned. A continued effort will be made 
toward understanding the chemical compo­ 
sition of the debris particles and manipu­ 
lation of the wettability of the particle 
surfaces to affect nucleation. 


Additional work is needed to more 
realistically describe rainout. 
Among 
the problems to be studied are what frac­ 


tion of debris particles can serve as con­ 
densation nuclei in a given precipitation 
system and how these particles compete 
with the natural condensation nuclei. 
If a 
large number of small particles result 
from a detonation, they would attach 
themselves to natural aerosol particles 
by virtue of their Brownian motion, and 
the removal of the contaminated natural 
aerosol particles must then be considered. 


The calculation of dry deposition can 
be improved by including the effect of 
wind shear on the traveling debris system. 
Another problem is the efficiency by 
which these particles attach to the earth's 
surface. 
Future work into this problem 
area is planned as well as an assessment 
of the global effects of wet and dry 
deposition. 


Fohl and Zalay 
suggested that the 
height of rise of two buoyant fireballs 
released simultaneously and in close 
proximity will be decreased by their 
mutual interaction. 
Further and more 
fundamental work is needed to evaluate 
and predict this phenomenon and to assess 
its implications on the scavenging of 
debris. 


In order to more effectively assess the 
consequences of the rainout problem, a 
climatological study for areas of interest 
is planned. Such information as the depth 
of convection cell size, the fraction of 
time that precipitating systems are pres­ 
ent, and the type of flow above precipitat­ 
ing systems coupled with improved rain­ 
out estimates, can be useful for military 
planning. 
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